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(57) Abstract 

The invention relates to a method for reducing noise and artefacts in digital subtraction imaging processes. A subtraction image data 
set is produced from two volume data sets, each volume data set consisting of a number of elements (voxels). A uni- or multidimensional 
vector (grey-scale value) is allocated to each voxel. A first m-dimensional limit vector, m corresponding to the number of dimensions of 
the vectors of the subtraction image data set, and a first proximity relation between the voxels of the subtraction image data set are defined. 
At least the n-te (n e N+) component of the m components of all of the vectors of the subtraction data set in each case is compared with 
the nth component of the first limit vector and when a predetermined first selection condition is fulfilled, are set at a predetermined first 
unit value. Each nth component of a voxel-allocated vector which does not fulfil the first selection condition is compared with the nth 
components of the vectors allocated to the voxels deemed closest to the voxel concerned according to the first proximity relation and when 
a predetermined second selection condition is fulfilled, is also set at the predetermined first unit value. 



(57) Zusammenfassung 

Bei einem Verfahren zur Rausch- und Artefaktreduktion bei digitalen Subtraktionsbildverfahren, bei welchen ein Subtraktionsbild- 
datensatz aus zwei Volumendatensatzen erzeugt wird, wobei jeder Volumendatensatz aus einer Vielzahl von Elementen (Voxeln) besteht 
und jedem Voxel ein ein- oder mehrdimensionaler Vektor (Grauwert) zugeordnet ist, wird ein erster m-dimensionaler Grenzvektor, wobei 
m der Anzahl der Dimensionen der Vektoren des Subtraktionsbilddatensatzes entspricht, und eine erste Nachbarschaftsbeziehung zwischen 
den Voxeln des Subtraktionsbilddatensatzes definiert, wenigstens die n-te (neN+) Komponente der jeweils m Komponenten aller Vektoren 
des Subtraktionsbilddatensatzes mit der n-ten Komponente des ersten Grenzvektors verglichen und bei Erfullung einer vorbestimmten ersten 
Auswahlbedingung auf einen vorbestimmten ersten Einheitswert gesetzt, wobei jede nicht die erste Auswahlbedingung erfiillende n-te Kom- 
ponente eines einem Voxel zugeordneten Vektors mit den n-ten Komponenten der dem Voxel gemaB der ersten Nachbarschaftsbeziehung 
benachbart liegenden Voxeln zugeordneten Vektoren verglichen und bei Erfullung einer vorbestimmten zweiten Auswahlbedingung eben- 
falls auf den vorbestimmten ersten Einheitswert gesetzt wird. 
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Verfahren zur Rausch- und Artefaktreduktion 
bei digitalen Subtraktionsbildverfahren 

Die Erfindung betrifft ein Verfahren zur Rausch- und Artefaktreduktion bei 
5 digitalen Subtraktionsbildverfahren, insbesondere bei digitalen 
Subtraktionsangiographieverfahren, bei welchen ein 
Subtraktionsbilddatensatz aus zwei Volumendatensatzen erzeugt wird, 
wobei jeder Volumendatensatz aus einer Vielzahl von Elementen (Voxeln) 
besteht und jedem Voxel ein ein- oder mehrdimensionaler Vektor 
10 (Grauwert) zugeordnet ist. 

Solche Verfahren sind bekannt, zum Beispiel aus dera Aufsatz vom K.-H. 
Englmeier et al. : "ein Verfahren zur computergestutzten, dreidimensionalen 
Subtraktionsangiographie mit Hilfe der Spiralrontgencomputertomographie" 
in: Radiologe (1996), 36: 360-364 oder dem Aufsatz von W. Lumboldt et 

15 al.: "CT-Subtraktionsangiographie - Ergebnisse mit einem automatisierten 
"elastischen" Subtraktionsalghorythmus", in: Radiologe (1997), 37 : 89-93. 
Sie dienen dazu, bei Subtraktionsbildverfahren, bei denen zwei 
Volumendatensatze voneinander "substrahiert" werden, Rauschen und 
Artefakte weitestmoglichst zu reduzieren. Die Ausgangsvolumendatensatze 

20 stammen dabei von einer elektromechanischen oder elektronischen 
Abtasteinrichtung, einem bildgebenden medizinischen Gerat, zum Beispiel 
einem Ultraschall- oder Rontgengerat, einem Computer-, Magnetresonanz- 
oder Positronenemissionstomographen oder wurden sonstwie erzeugt, zum 
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Beispiel aus gemessenen Daten berechnet. Beide Volumendatensatze 
enthalten Informationen iiber ein- und dasselbe Aufnahmegebiet, 
unterscheiden sich jedoch zum Beispiel hinsichtlich bestimmter 
Aufnahmeparameter, zum Beispiel der Leistung, mit der ein Rontgengerat 
5 bei der Aufnahme der beiden Volumendatensatze betrieben wurde. Sie 
konnen sich aber auch dadurch unterscheiden, daft physikalische und/oder 
chemische Eigenschaften des sich in dem Aufnahmegebiet befindenden, zu 
untersuchenden Korpers verandert wurden, wie zum Beispiel bei der 
Subtraktionsangiographie, bei der mit derselben Rontgenleistung zwei 
10 Volumendatensatze ein und desselben Korperteils, insbesondere des 
menschlichen Kopfes, erzeugt werden, wobei die Aufnahme jedoch einmal 
mit und einmal ohne vorherige Verabreichung eines Kontrastmittels 
gemacht wird. 

Schliefilich konnen die beiden Volumendatensatze auch von ganzlich 
15 unterschiedlichen Aufnahmeeinrichtungen stammen, zum Beispiel vom 
einem Computertomographen, der einen guten Knochen-Weichteil-Kontrast 
liefert, und einem Magnetresonanztomographen, der einen guten Weichteil- 
Weichteil-Kontrast liefert. Durch die Subtraktion der beiden 
Volumendatensatze erhofft man sich eine Verbesserung der in den 
20 Volumendatensatzen enthaltenen Informationen. So konnen zum Beispiel 
bei der Subtraktionsangiographie, bei der es um die Sichtbarmachung von 
GefaBen geht, die hier nur storenden Informationen iiber vorhandene 
Knochen weitestgehend eliminiert werden. 
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Bei den Ausgangs-Volumendatensatzen kann es sich um in einem 
Aufnahmeschritt aufgenommene Volumendatensatze handeln, wie sie zum 
Beispiel von bestimmten medizinischen bildgebenden Geraten, zum Beispiel 
3D-Ultraschallwandlern erzeugt werden, oder um aus einer Anzahl von 
5 zweidimensionalen Schnittbildern, wie sie zum Beispiel von den meisten 
Computertomographen oder iiblichen 2D-Ultraschallwandlern geliefert 
werden, zusammengesetzte Volumendatensatze handeln. Die beiden 
Ausgangsvolumendatensatze besitzen in der Regel, insbesondere nariirlich 
immer dann, wenn sie von ein und demselben Aufhahmegerat stammen, die 

10 gleichen Dimensioned bestehen also zum Beispiel aus 136 x 256 x 256 
Voxeln. Es konnen aber auch Volumendatensatze zueinander in Beziehung 
gesetztund voneinander subtrahiert werden, die nicht aus derselben Anzahl 
von Voxeln bestehen, was in der Regel dann auftritt, wenn die 
Volumendatensatze nicht von einem einzigen Aufhahmegerat aufgenommen 

15 wurden. Dabei wird der Aufnahmeraum in an sich bekannter Weise in 
kleine Raumeinheiten unterteilt und die in eine bestimmte Raumeinheit 
fallenden Voxel des einen Volumendatensatzes werden nach entsprechender 
Interpolation von den in dieselbe Raumeinheit fallenden Voxeln des 
anderen Volumendatensatzes "subtrahiert" , wobei natiirlich nicht Voxel, die 

20 ja nur den Ort beschreiben, an dem eine bestimmte Information 
aufgenommen wurde, sondern die jedem Voxel in Form eines ein- oder 
mehrdimensionalen Vektors zugeordneten Werte voneinander subtrahiert 
werden. 

Den einzelnen Voxeln des Volumendatensatzes sind also Vektoren 
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zugeordnet, die im einfachsten Fall, wenn es sich urn Grauwertbilder 
handelt, eindimensional (skalar) sind und deren Komponenten in der Regel 
nur aus einer ganzen Zahl aus einem vorbestimmten Wertebereich, zum 
Beispiel [0,255], mit einer vorbestimmten Anzahl, zum Beispiel 256 
5 moglicher Zahlenwerte besteht, die dann moglichen Graustufen des 
jeweiligen Volumendatensatzes entsprechen, wobei zum Beispiel der 
niedrigste Wert fur ein weifies, der hochste Wert fur ein schwarzes Voxel 
steht, wobei aber auch beliebige andere Zuordnungen moglich sind. 

Den Voxeln konnen anstelle eindimensionaler auch, wie zum Beispiel bei 
10 Farbbildern, dreidimenionale Vektoren zugeordnet sein, deren 
Komponenten den moglichen Intensitatsniveaus dreier Grundfarben, zum 
Beispiel den Farben rot, griin, blau oder zyan, magenta und gelb 
entsprechen. 

SchlieBlich konnen den Voxeln auch Vektoren mit mehr ais drei 
15 Dimensionen zugeordnet sein, zum Beispiel 20-dimensionale wie bei 
einigen Radarbeobachtungssatelliten, deren 20 verschiedenen Komponenten 
dann zum Beispiel den moglichen Intensitatsniveaus verschiedener 
Meligroflen entsprechen. 

Die in den Volumendatensatzen und insbesondere dem daraus erzeugten 
20 Subtraktionsbilddatensatz enthaltenen Informationen konnen in an sich 
bekannter Weise angezeigt werden, zum Beispiel in Form einer plastischen 
Pseudo-3D-Darstellung auf einem Bildschirm oder einem Drucker 
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ausgegeben werden, wozu die unterschiedlichsten Verfahren, zum Beispiel 
sogenannte Raytracing-Verfahren oder die sogenannte maximale 
Intensitatsprojektion verwendet werden. 

Nun sind nicht nur die Ausgangs-Volumendatensatze rait Rauschen und je 
5 nach Aufnahmeverfahren - auch mit Artefakten behaftet, auch beim 
"Subtrahieren" treten aus unterschiedlichsten Griinden neues Rauschen und 
neue Artefakte auf. Insbesondere bei der Subtraktionsangiographie beruht 
ein Teil der auftretenden Artefakte auf translatorischen Verschiebungen des 
untersuchten Korpers zwischen der Aufnahme der beiden 

10 Volumendatensatze. Das in dem genannten Aufsatz von Englmeier 
beschriebene Verfahren ermittelt solche Translationen unter Verwendung 
der sogenannten "Cepstrumfiltermethode". Es kann zwar einen Teil der 
auftretenden Artefakte beseitigen, ist jedoch aufgrund der bei diesem 
Verfahren durchzufiihrenden Interpolationen zum einem 

15 rechenzeitaufwendig und zum anderen mit einem Informationsverlust 
verbunden. Zudem erscheint fraglich, ob das Verfahren Artefakte, die 
aufgrund einer fehlenden Synchronitat von Rontgenrohre und Rontgentisch 
auftreten, iiberhaupt beseitigen kann. Das in dem genannten Aufsatz von 
Lumboldt beschriebene Verfahren kann dem gegeniiber sowohl 

20 Translationen als auch Verzerrungen korrigieren, wozu in den 
Volumendatensatzen kongruente Voluminaflachen, Kanten und Punkte 
elastisch zur Deckung gebracht werden. Dieses sogenannte "elastische 
Matchen" ist jedoch ebenfalls sehr rechenzeitaufwendig, und ist nicht nur 
mit einem Informationsverlust aufgrund der auch bei diesem Verfahren 
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auszuftihrenden Interpolationen, sondern auch mit einem Verlust der gerade 
bei der Subtraktionsangiographie sehr hohen geometrischen Exaktheit 
verbunden. 

Davon ausgehend liegt der Erfmdung die Aufgabe zugrunde, ein Verfahren 
5 zur Rausch- und Artefaktreduktion bei Subtraktionsbildverfahren 
anzugeben, das zum einen Rauschen und Artefakte weitestgehend ohne 
Informationsverlust reduziert und das zum anderen eine moglichst geringe 
Komplexitat aufweist, so dafi das Verfahren sehr schnell arbeitet. 

Die Aufgabe wird gelost von einem Verfahren der eingangs genannten Art, 
10 bei welchem ein erster m-dimensionaler Grenzvektor, wobei m der Anzahl 
der Dimensionen der Vektoren des Subtraktionsbilddatensatzes entspricht, 
und eine erste Nachbarschaftsbeziehung zwischen den Voxeln des 
Subtraktionsbilddatensatzes definiert werden, wenigstens die n-te (wobei n 
Element der positiv naturlichen Zahlen ist) Komponente der jeweils m- 
15 Komponenten aller Vektoren des Subtraktionsbilddatensatzes mit der n-ten 
Komponente des ersten Grenzvektors verglichen und bei Erfullung einer 
vorbestimmten ersten Auswahlbedingung auf einen vorbestimmten ersten 
Einheitswert gesetzt wird und wobei jede nicht die erste Auswahlbedingung 
erfullende n-te Komponete eines einem Voxel zugeordneten Vektors mit 
20 den n-ten Komponenten der dem Voxel gemaii der ersten 
Nachbarschaftsbeziehung benachbart liegenden Voxeln zugeordneten 
Vektoren und bei Erfullung einer vorbestimmten zweiten 
Auswahlbedingung ebenfalls auf den vorbestimmten ersten Einheitswert 
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gesetzt wird. Dabei enthalt die erste Auswahlbedingung vorzugsweise eine 
Abfrage, ob der Wert der jeweils mit der entsprechenden Komponente des 
ersten oder zweiten Grenzvektors verglichenen Komponente des jeweiligen 
Vektors iiber- oder unterhalb des Wertes der entsprechenden Komponente 
5 des jeweiligen Grenzvektors liegt. 

Die zweite Auswahlbedingung kann vorteilhaft eine Ahnlichkeitsabfrage 
enthalten, ob einem oder mehreren der gemaft der ersten 
Nachbarschaftsbeziehung dem jeweils betrachteten Voxel benachbart 
liegenden Voxeln Vektoren zugeordnet sind, deren Komponenten der 
10 jeweils betrachteten Komponente des dem betrachteten Voxel zugeordneten 
Vektors ahnlich sind. 

Die Erfindung beruht also auf dem Grundgedanken, daft nach der 
Subtraktion die verbliebenen Artefakte zufallig verteilt sind und keine 
grofieren zusammenhangenden Strukturen bilden. Stammen beide 

15 Datensatze von ein und demselben Aufnahmegerat, ist diese Annahme 
jedenfalls immer dann gerechtfertigt, wenn sich das Objekt zwischen den 
Aufhahmen der beiden Datensatze gar nicht oder nur sehr wenig bewegt 
hat. Insbesondere bei medizinischen Bilddatensatzen wird also in der Regel 
auf eine ausreichende Immobilisierung in an sich bekannter Weise zu 

20 achten sein. Trifft die genannte Annahme zu, konnen alle keine 
zusammenhangenden Strukturen bildenden Voxel oder Ansammlungen von 
Voxeln (Cluster) auf einen vorbestimmten ersten Einheitswert, zum 
Beispiel Null gesetzt werden und erscheinen bei einer anschlieftenden 
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Darstellung des Subtraktionsbilddatensatzes zum Beispiel mittels der 
erwahnten maximalen Intensitatsprojektion nicht und sind gewissermafien 
"durchsichtig". Es sei an dieser Stelle betont, daB wenn hier davon 
gesprochen wird, da6 "ein Voxel auf einen bestimmten Einheitswert" 
5 gesetzt wird, damit natiirlich gemeint ist, dafl die jeweils betrachtete 
Komponente des dem Voxel zugeordnetem Vektors auf diesen Einheitswert 
gesetzt wird. Bei reinen Grauwertbildern sind die den Voxeln zugeordneten 
Vektoren natiirlich eindimensional und jedem Voxel ist nur ein einziger, 
skalarer Wert zugeordnet. Aus historischen Griinden werden allerdings 

10 auch mehrdimensionale Vektoren oftmals als "Grauwert" der Voxel 
bezeichnet, wobei das erfindungsgemaBe Verfahren bei 
Volumendatensatzen aus Voxeln, denen mehrdimensionale Vektoren 
zugeordnet sind, dann zum einen so durchgefuhrt werden kann, daB jeweils 
nur eine bestimmte Komponente der Vektoren, die zum Beispiel dem 

15 Meliwert fur das Absorptionsverhalten in dem Voxel fur Strahlung einer 
bestimmten Intensitat entspricht, betrachtet wird. Das erfindungsgemafie 
Verfahren kann aber auch so durchgefuhrt werden, daJi nacheinander oder 
abwechselnd mehrere oder alle Komponenten der Vektoren daraufhin 
untersucht werden, ob sie Rauschen und/oder Artefakte reprasentieren, um 

20 dann auf den bzw. die vorbestimmten Einheitswerte (in der Regel Null) 
gesetzt zu werden. 



Eine besonders einfache Durchfuhrungsform des erfindungsgemafien 
Verfahrens lafit sich vorteilhaft dann realisieren, wenn die Voxel der 
beiden Volumendatensatze und der Schnittbildsatze jeweils einem Pixel 
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eines aus einer Vielzahl von Pixeln bestehenden Schnittbildes einer den 
Volumen- und den Bilddatensatz bildenden Anzahl von Schnittbildern 
entsprechen, was zum Beispiel bei den meisten bildgebenden Verfahren in 
der Medizin, insbesondere der Subtraktionsangiographie der Fall ist. Dann 
5 kann namlich die erste Nachbarschaftsbeziehung derart definiert werden, 
dafi zu jedem Pixel einer Schicht nur ein Pixel aus anderen, an die Schicht 
des jeweiligen Pixels angrenzenden Schichten, nicht aber Pixel aus 
derselben Schicht als benachbart angesehen werden. Die Erzeugung des 
Subtraktionsbilddatensatzes wird namlich bei solchen aus einer Anzahl von 

10 Schnittbildern bestehenden Volumendatensatzen der einfachheithalber so 
vorgenommen, daft jeweils einander raumlich entsprechende Schichten aus 
beiden Volumendatensatzen voneinander subtrahiert werden, so daft 
Subtraktionsartefakte in jeder Schicht gesondert entstehen und in der Regel 
keine sich fiber mehrere Schichten hinwegziehenden zusammenhangenden 

15 Cluster bilden. Man kann daher diesen Schritt des erfindungsgemafien 
Verfahrens als "Kontinuitatscheck" bezeichnen, bei dem also gepriift wird, 
ob ein aufgrund seiner zum Beispiel oberhalb bestimmter, in dem ersten 
Grenzvektor festgelegter Grenzwerte liegenden "Grauwerte" als 
Artefaktvoxel in Betracht kommendes Voxel Entsprechungen in seiner 

20 Nachbarschaft bzw. benachbarten Schichten hat. Hat es diese 
Entsprechungen -hat also der Grauwert iiber mehrere Voxel hinweg eine 
gewisse Kontinuitat - so wird davon ausgegangen, dafi es sich nicht um 
einen zufallig entstandenen Artefakt, sondernum tatsachlich interessierende 
Informationen handelt und die Komponenten der den entsprechenden 

25 Voxeln zugeordneten Vektoren bleiben unverandert. 
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In vorteilhafter Weiterbildung der Erfindung kann wenigstens ein zweiter 
m-dimensionaler Grenzvektor definiert werden und vor der Erzeugung des 
Subtraktionsbilddatensatzes kann wenigstens eine Komponente der Vektoren 
wenigstens eines der beiden Volumendatensatze mit der entsprechenden 

5 Komponente des Grenzvektors verglichen und bei Erfullung einer dritten 
Auswahlbedingung auf einen vorbestimmten, nicht notwenigerweise mit 
dem ersten Einheitswert identischen zweiten Einheitswert gesetzt werden. 
Dies entspricht dem, was man in der Bildverarbeitung allgemein als 
"Glattung" bezeichnet und dient dazu, bestimmte, unter oder oberhalb eines 

10 gewissen Schwellenwertes liegende "Grauwerte" zu unterdriicken, so dali 
bei der anschliefienden Subtraktion der beiden Volumendatensatze 
tatsachlich nur diejenigen Voxel voneinander "subtrahiert" werden, die 
uberhaupt dafur in Frage kommen, relevante Informationen zu enthalten. 

Bei einer bevorzugten Ausfuhrungsform der Erfindung ist vorgesehen, dali 
15 bei der Erzeugung des Subtraktionsbilddatensatzes nur diejenigen Voxel der 
Volumendatensatze beriicksichtigt werden, deren zugeordnete Vektoren 
Komponenten besitzen, die in vordefinierte erste, nicht notwendigerweise 
fur beide Volumendatensatze identische m-dimensionale Wertebereiche 
fallen. In der Regel ist es dabei vorteilhaft vorzusehen, dali die 
20 Wertebereiche fur jeden der Volumendatensatze verschieden sind. Man 
kann diesen Verfahrensschritt als "parametrisierte Subtraktion" bezeichnen 
und sie dient ahnlich wie die zuvor beschriebene Glattung dazu, nur 
diejenigen Voxel bei der Erzeugung des Subtraktionsbilddatensatzes zu 
beriicksichtigen, die tatsachlich relevante Informationen enthalten, was 
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nicht nur die Bildqualitat erheblich verbessert, sondern was auch die 
Komplexitat der auszufuhrenden Subtraktionen und - je nach Umsetzung - 
auch die Anzahl der auszufuhrenden Subtraktionen erheblich verringertund 
es damit ermoglicht, das erfindungsgemaBe Verfahren besonders schnell 
5 auszufiihren. 

Je nach Art des Subtraktionsbildverfahrens, bei dem das erfindungsgemalie 
Verfahren angewandt werden soil, kann es vorkommen, daft bei der zuvor 
beschriebenen "parametrisierten Subtraktion" noch Voxel beriicksichtigt 
wurden und zu entsprechenden Artefakten fuhren, die an sich nicht weiter 

10 inter essieren. Diese Voxel konnen in einem sogenannten 
"Ausdunnungsschritt" nach der Erzeugung des Subtraktionsbilddatensatzes 
und vor dem Vergleich der Vektoren des Subtraktionsbilddatensatzes mit 
dem ersten Grenzvektor dadurch ausgesondert werden, dafi eine nicht 
notwendigerweise mit der ersten Nachbarschaftsbeziehung identische zweite 

15 Nachbarschaftsbeziehung definiert wird und daiJ die Voxel des 
Subtraktionsbilddatensatzes, denen Vektoren mit Komponenten innerhalb 
eines vorbestimmten zweiten m-dimensionalen Wertebereichs zugeordnet 
sind, ausgewahlt werden und dafl wenigstens eine Komponente des Vektors 
jedes ausgewahlten Voxels dann auf einen vorbestimmten dritten 

20 Einheitswert (in der Regel Null) gesetzt wird, wenn bei der Erzeugung 
dieser ausgewahlten Voxeln beriicksichtigten Voxeln eines der beiden 
Volumendatensatze gemafl der zweiten Nachbarschaftsbeziehung Voxel 
benachbart liegen, die bei der Erzeugung des Subtraktionsbilddatensatzes 
aufgrund entweder des Uber- oder des Unterschreitens des zugeordneten 
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Wertebereichs nicht beriicksichtigt wurden. Dieser vorteilhafte 
"Ausdunnungsschritt" des erfindungsgemafien Verfahrens wird - ebenso 
wie der "Kontinuitatscheck" und die "parametrisierte Subtraktion" - 
nachfolgend am Beispiel der Anwendung des erfindungsgemafien 
5 Verfahrens bei der Subtraktionsangiographie unter Bezugnahme auf die 
beigefugten Figuren naher erlautert. Es zeigen, 

Figuren 1 und 2 zwei von einem Computertomographen jeweils an der 
selben Stelle aufgenommene Schnittbilder eines 
menschlichen Schadels, jedoch einmal ohne (Figur 2) 
10 und einmal mit (Figur 1) vorheriger Verabreichung 



eines Kontrastmittels, 



Figur 3 



das Ergebnisbild einer reinen "voxelweisen" Subtraktion 
des in Figur 2 gezeigten Schnittbildes von dem in Figur 
1 gezeigten Schnittbild, 



15 



Figur 4 



das Ergebnis einer nach dem erfindungsgemafien 
Verfahren ausgefuhrten "parametrisierten Subtraktion" 
der Schnittbilder gemafi Figur 1 und Figur 2 
voneinander, 



20 



Figur 5 



eine plastisch wirkende maximale Intensitatsprojektion 
eines kompletten, mittels der parametrisierten 
Subtraktion erzeugten Subtraktionsbilddatensatzes, 
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Figur 6 das Schnittbild gemafi Figur 4, jedoch nach dera 

Ausfuhren des Ausdiinnungsschritts und des 
Kontinuitatschecks und 

Figur 7 eine maximale Intensitatsprojektion des kompletten, aus 

einer Anzahl von Schnittbildern wie den Figur 6 
gezeigten bestehenden Subtraktionsbilddatensatzes. 

In den Figuren 1 und 2 sind zwei von einem Computertomographen jeweils 
am selben Ort aufgenommene Schnittbilder eines menschlichen Schadels 
gezeigt, wobei jedoch vor der Aufhahme des in Figur 1 gezeigten 
Schnittbildes dem untersuchten Patienten ein Kontrastmittel verabreicht 
wurde, so dafi die durchbluteten Gefafie in Figur 1 heller erscheinen als in 
Figur 2. Im Abstand weniger Millimeter aufgenommene Schnittbilder der 
in Figuren 1 und 2 gezeigten Art bilden die Ausgangs-Volumendatensatze, 
aus denen nun Subtraktionsbilddatensatze erzeugt werden sollen. Dabei 
besteht ein solcher Voluraendatensatz zum Beispiel aus 136 Schnittbildern, 
die jeweils aus zum Beispiel 256 x 256 Pixeln bestehen. Jedem Pixel (wenn 
man ein einziges Schnittbild betrachtet) bzw. jedem Voxel (wenn man den 
gesamten Volumendatensatz betrachtet) ist ein im Falle der 
Computertomographie eindimensionaler Vektor zugeordnet, dessen einzige 
Komponente nur ganzzahlige Werte aus einem vorbestimmten 
Wertebereich, zum Beispiel [0,255] annehmen kann, wobei der 
angenommene Wert dem Absorptionsverhalten fur Rontgenstrahlung einer 
bestimmten Wellenlange, namlich der Wellenlange des verwendeten 
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Computertomographen entspricht. 

Bei alien hier gezeigten Figuren ist dabei die Umsetzung der gemessenen 
Absorptionswerte in bildlich darstellbare Grauwerte so gewahlt, daB 
Bildelemente, in denen eine starkere Absorption der Strahlung gemessen 
5 wurde, heller erscheinen als Bildelemente, die die Strahlung nur wenig 
oder gar nicht absorbieren. Die Schadeldecke des Patienten erscheint daher 
in den Schnittbildern der Figuren 1 und 2 als dicke, ringformige, weiBe 
Linie. 

Da bei Computertomographen der Weichteil-Weichteil-Kontrast bekannter- 
10 weise schlecht ist, wird bei der Subtraktionsangiographie ein Kontrastmittel 
verabreicht, so daB die Gefafle bei der Aufnahme mit Kontrastmittel heller 
erscheinen (Figur 1) als bei der Aufnahme derselben GefaBe ohne 
Kontrastmittel (Figur 2). Da die Helligkeit der hier nicht interessierenden 
Knochen bei beiden Aufhahmen an sich identisch sein imiBte, da sie nicht 
15 durch das Kontrastmittel beeinflufit wird, sollten die Knochen und weiteren 
nicht interessierenden Gebiete bei der Subtraktion des an der selben Stelle 
aufgenommenen Schnittbildes gemaB Figur 2 von dem Schnittbild gemafi 
Figur 1 wegfallen und es sollten in dem Differenzbild nur die durchbluteten 
und damit durch das Kontrastmittel hervorgehobenen GefaBe 
20 zuriickbleiben. 

Wie Figur 3 zeigt, ist dies nicht der Fall und man erkennt neben dem 
Rauschen noch die typischen Knochenartefakte. Diese stammen nicht von 
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einer Translation, da sich sonst sowohl an der frontalen als auch an der 
okzipitalen Kalottenseite die typischen hellen und dunklen Streifen zeigten. 

Wie in Figur 4 gezeigt, kann die erfmdungsgemafl ausgefiihrte 
"parametrisierte Subtraktion" diese Artefakte bereits deutlich verringern. 

5 Bei der "parametrisierten Subtraktion" wird - angewendet auf die 
Subtraktionsangiographie - davon ausgegangen, da/5 die interessierenden 
Gefafie in einem ohne Kontrastmittel aufgenommenen Schnittbild nur 
Absorptionswerte innerhalb eines bestimmten ersten Wertebereichs besitzen 
konnen. Es wird also zum Beispiel angenommen, dafi alle helleren Gebiete 

10 in dem Schnittbild von Knochen stammen. In analoger Weise wird 
angenommen, dafl die interessierenden Gefafle nach Verabreichung des 
Kontrastmittels in einem Schnittbild nur Werte innerhalb eines bestimmten 
zweiten Wertebereichs annehmen konnen, der bei diesem Beispiel natiirlich 
oberhalb des ersten Wertebereichs liegen muli, da das Kontrastmittel gerade 

15 dazu fiihrt, dafl die Gefafle heller erscheinen. Bei der parametrisierten 
Subtraktion werden nun iiberhaupt nur diejenigen Voxel beriicksichtigt, 
denen Werte zugeordnet sind, die innerhalb des fur den jeweiligen 
Volumendatensatz zugeordneten Wertebereichs liegen. Wie Figur 4 zeigt, 
sind die Artefakte bereits signifikant reduziert, jedoch noch nicht 

20 zufriedenstellend unterdriickt. Eine maximale Intensitatsproj ektion eines aus 
einer Anzahl von Schnittbildern gemaB Figur 4 bestehenden 
Subtraktionsbilddatensatzes (Figur 5) zeigt noch Storungen. Werden nun 
der Ausdiinnungsschritt und der Kontinuitatscheck durchgefuhrt, so erhalt 
man diagnostisch wertvolle Schnittbilder, wie in Figur 6 gezeigt, wobei ein 
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ganzer Satz solcher Schnittbilder dann zum Beispiel mit der genannten 
maximalen Intensitatsprojektion als dreidimensionales Bild dargestellt 
werden kann (Figur 7). Dabei beruht der sogenannte Ausdunnungsschritt 
auf der Annahme, dafi speziell bei der Computertomographie Knochen in 
5 ihren Randbereichen aufgrund von Partialvolumeneffekten 
Absorptionswerte besitzen konnen, die deutlich geringer sind, als im 
Zentrum des Knochens, so dafi die Randbereiche der Knochen zu Voxeln 
fuhren, denen Grauwerte zugeordnet sind, die innerhalb des ersten oder 
zweiten Wertebereichs liegen und die deshalb bei der parametrisieren 

10 Subtraktion nicht erfafit werden konnen, da ihre Helligkeit der Helligkeit 
der interessierenden Gefafie entspricht. Ein Nichtberiicksichtigen dieser 
Voxeln bei der parametrisierten Subtraktion konnte also nur dadurch 
erreicht werden, dafi die Wertebereiche entsprechend verkleinert wurden, 
was dann aber zur Konsequenz hatte, dafi auch die eigentlich 

15 interessierenden Gefafie herausfielen. Beim sogenannten 
Ausdunnungsschritt wird nun iiberpriift, ob in der Nachbarschaft solcher 
Voxel, die bei der parametrisierten Subtraktion beriicksichtigt wurden, 
eventuell Voxel liegen, die bei der Subtraktion aufgrund zu hoher 
Grauwerte als Knochen identifiziert und bei der Subtraktion nicht 

20 beriicksichtigt wurden. 1st dies der Fall, so wird davon ausgegangen, dafi 
auch die zu diesen nicht beriicksichtigten Voxeln benachbart liegenden, 
jedoch nicht ganz so hell erscheinenden Voxel Teile des nicht 
interessierenden Knochens, namlich dessen stets nicht so stark 
absorbierende Randbereiche reprasentieren. Die auf diese Weise als zum 

25 Knochen gehorig identifizierten Voxel werden dann entsprechend behandelt 
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und zum Beispiel, wie in den Figuren 6 und 7 gezeigt, auf "Null" gesetzt, 
so dali sie in den Bildern nicht mehr in Erscheinung treten und die reinen, 
eigentlich interessierenden GefaBstrukturen deutlich sichtbar werden. Es sei 
an dieser Stelle betont, daB dieser Ausdiinnungsschritt natiirlich theoretisch 
5 auch vor der Subtraktion ausgefuhrt werden konnte, so daB bereits in den 
urspriinglichen Volumendatensatzen diejenigen Voxel identifiziert werden 
konnten, die Knochenrandbereiche reprasentieren. Da jedoch in den 
urspriinglichen Datensatzen entsprechend viele "Knochenvoxel" enthalten 
sind, ware eine solche Vorgehensweise sehr aufwendig. Besser ist es daher, 
10 zunachst durch die einfachere Subtraktionsoperation bereits einen GroBteil 
der Knochenvoxel zu entfernen, so daB im Ausdiinnungsschritt nur noch 
relativ wenige mogliche Knochenrandbereichsvoxel identifiziert und 
ausgedunnt werden mussen, was vorteilhaft eine besonders schnelle 
Ausfuhrung des Verfahrens sicherstellt. 
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Patentanspriiche: 

1 . Verfahren zur Rausch- und Artefaktreduktion bei digitalen Subtrak- 
tionsbildverfahren, insbesondere bei digitalen Subtrak- 
tionsangiographieverfahren, bei welchen em Subtraktionsbild- 
datensatz aus zwei Volumendatensatzen erzeugt wird, wobei jeder 
Volumendatensatz aus einer Vielzahl von Elementen (Voxeln) besteht 
und jedem Voxel ein ein- oder mehrdimensionaler Vektor (Grauwert) 
zugeordnet ist, dadurch gekennzeichnet, dafi ein erster m- 
dimensionaler Grenzvektor, wobei m der Anzahl der Dimensionen 
der Vektoren des Subtraktionsbilddatensatzes entspricht, und eine 
erste Nachbarschaftsbeziehung zwischen den Voxeln des 
Subtraktionsbilddatensatzes definiert werden, daft wenigstens die n-te 
(n G N + ) Komponente der jeweils m Komponenten aller Vektoren 
des Subtraktionsbilddatensatzes mit der n-ten Komponente des ersten 
Grenzvektors verglichen und bei Erfullung einer vorbestimmten 
ersten Auswahlbedingung auf einen vorbestimmten ersten 
Einheitswert gesetzt wird, daft jede nicht die erste Auswahlbedingung 
erfullende n-te Komponente eines einem Voxel zugeordneten Vektors 
mit den n-ten Komponenten der dem Voxel gemaft der ersten 
Nachbarschaftsbeziehung benachbart liegenden Voxeln zugeordneten 
Vektoren verglichen und bei Erfullung einer vorbestimmten zweiten 
Auswahlbedingung ebenfalls auf den vorbestimmten ersten 
Einheitswert gesetzt wird. 
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Verfahren nach Anspruch 1, dadurch gekennzeichnet, daft, 
wenigstens ein zweiter m-dimensionaler Grenzvektor definiert wird 
und dafi vor der Erzeugung des Subtraktionsbilddatensatzes 
wenigstens eine Komponente der Vektoren wenigstens eines der 
beiden Volumendatensatze mit der entsprechenden Komponente des 
Grenzvektors verglichen und bei Erfullung einer dritten Auswahl- 
bedingung auf einen vorbestimmten, nicht notwendigerweise mit dem 
ersten Einheitswert identischen zweiten Einheitswert gesetzt wird. 

Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, dafi 
bei der Erzeugung des Subtraktionsbilddatensatzes nur diejenigen 
Voxel der Volumendatensatze berucksichtigt werden, deren 
zugeordnete Vektoren Komponenten besitzen, die in vordefinierte 
erste, nicht notwendigerweise fur beide Volumendatensatze iden- 
tische m-dimensionale Wertebereiche fallen. 

Verfahren nach Anspruch 3, dadurch gekennzeichnet, dafi die 
Wertebereiche fur jeden Volumendatensatz verschieden sind. 

Verfahren nach Anspruch 3 oder 4, dadurch gekennzeichnet, dafi 
eine nicht notwendigerweise mit der ersten Nachbarschaftsbeziehung 
identische zweite Nachbarschaftsbeziehung definiert wird, daB vor 
dem Vergleich der Vektoren des Subtraktionsbilddatensatzes mit dem 
ersten Grenzvektor die Voxel des Subtraktionsbilddatensatzes, denen 
Vektoren mit Komponenten innerhalb eines vorbestimmten zweiten 
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m-dimensionalen Wertebereichs zugeordnet sind, ausgewahlt werden 
und dafl wenigstens eine Komponente des Vektors jedes 
ausgewahlten Voxels dann auf einen vorbestimmten dritten 
Einheitswert gesetzt wird, wenn den bei der Erzeugung dieser 
ausgewahlten Voxeln berucksichtigten Voxeln eines der beiden 
Volumendatensatze gemafi der zweiten Nachbarschaftsbeziehung 
Voxel benachbart liegen, die bei der Erzeugung des Sub- 
traktionsbilddatensatzes aufgrund entweder des Uber- oder des 
Unterschreitens des zugeordneten Wertebereichs nicht beriicksichtigt 
wurden. 

6. Verfahren nach Anspruch 5, dadurch gekennzeichnet, dafi die 
Auswahl der Voxel des Subtraktionsbilddatensatzes, denen Vektoren 
mit Komponenten innerhalb eines vorbestimmten zweiten m- 
dimensionalen Wertebereichs zugeordnet sind, und das Setzen der 
wenigstens einen Komponente des Vektors jedes ausgewahlten 
Voxels auf einen vorbestimmten dritten Einheitswert, wenn den bei 
der Erzeugung dieser ausgewahlten Voxeln berucksichtigten Voxeln 
eines der beiden Volumendatensatze gemafi der zweiten 
Nachbarschaftsbeziehung Voxel benachbart liegen, die bei der 
Erzeugung des Subtraktionsbilddatensatzes aufgrund entweder des 
Uber- oder des Unterschreitens des zugeordneten Wertebereichs 
nicht beriicksichtigt wurden, nach der Erzeugung des 
Subtraktionsbilddatensatzes erfolgt. 
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7. Verfahren nach einem der vorangehenden Anspriiche, wobei die 

Position jedes Voxels durch drei Indizes i, j, k £ N+ eindeutig 
beschrieben werden kann, dadurch gekennzeich.net, dafi die 
erste und/oder die zweite Nachbarschaftsbeziehung eine vorbe- 
5 stimmte Anzahl von Inkrementationen und/oder 

Dekrementationen eines oder mehrerer Indizes enthalt. 

8. Verfahren nach einem der vorangehenden Anspriiche, dadurch 

gekennzeichnet, da/Sjedem Voxel ein eindimensionaler Vektor 
zugeordnet ist, dessen Komponente nur ganzzahlige Werte aus 
10 einem vorgegebenen Bereich mdglicher Werte annehmen kann. 



9. Verfahren nach Anspruch 8, dadurch gekennzeichnet, dafi das 
wenigstens einer, vorzugsweise alle Einheitswerte entweder 
dem niedrigsten oder dem hochsten Wert des Bereichs 
mdglicher Werte entsprechen. 

15 10. Verfahren nach einem der vorangehenden Anspriiche, dadurch 
gekennzeichnet, dafi die erste und/oder dritte Auswahlbe- 
dingung eine Abfrage enthalt, ob der Wert der jeweils mit der 
entsprechenden Komponente des ersten oder zweiten 
Grenzvektors verglichenen Komponente des jeweiligen 

20 Vektors iiber- oder unterhalb des Wertes der entsprechenden 

Komponente des jeweiligen Grenzvektors liegt. 
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11. Verfahren nach einem der vorangehenden Anspriiche, dadurch 
gekennzeichnet, dafi die zweite Auswahlbedingung eine Ahnlich- 
keitsabfrage enthalt, ob einem oder mehreren der gemafl der ersten 
Nachbarschaftsbeziehung dem jeweils betrachteten Voxel benachbart 
liegenden Voxeln Vektoren zugeordnet sind, deren Komponenten der 
jeweils betrachteten Komponente des dem betrachteten Voxel 
zugeordneten Vektors ahnlich sind. 

12. Verfahren nach einem der vorangehenden Anspriiche, wobei jedes 
Voxel der beiden Volumendatensatze und des Subtraktionsbild- 
datensatzes jeweils einem Pixel eines aus einer Vielzahl von Pixeln 
bestehenden Schnittbildes einer die Volumendatensatze und den 
Subtraktionsbilddatensatz bildenden Anzahl von Schnittbildern 
entspricht, dadurch gekennzeichnet, dafi die erste Nachbar- 
schaftsbeziehung derart definiert ist, daft zu jedem Pixel einer 
Schicht nur Pixel aus anderen, angrenzenden Schichten, nicht aber 
Pixel aus derselben Schicht als benachbart angesehen werden. 
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